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Introduction 

One of the most reported problem related to the CMNS (condensed matter nuclear science, or 
LENR), is the low probability of Coulomb barrier tunneling. It is supposed by standard physics 
that tunneling is only possible at high enough energy (by solving Gamow function). 

However, a recent study by A. Takahashi (2008, 2009) and experiment by Arata etc. (2008) 
seem to suggest that it is not impossible to achieve a working experiment to create the CMNS 
process. 

In accordance with Takahashi’ s EQPET/TSC model [1][2][3], the proposed study will find out 
some analytical and numerical solutions to the problem of barrier tunneling for cluster 
deuterium, in particular using Langevin method to solve the time-independent Schrodinger 
equation. It is hoped that the result can answer some of these mysteries. 

One of the results of recent experiments is the lack of signature of D-D reaction as in standard 
fusion process; this is part of the reason to suggest that D-D fusion doesn’t take place [1], 
However, Takahashi suggests new possible reaction in the context of cluster deuterium, called 
4D fusion [1][2][3], this mechanism seems to enable reaction at low temperature (CMNS). His 
result (2009) can be summarized as follows: 

“The ultimate condensation is possible only when the double Platonic symmetry of 
4D/TSC is kept in its dynamic motion. The sufficient increase (super screening) of barrier 
factor is also only possible as far as the Platonic symmetric 4D/TSC system is kept. 
Therefore, there should be always 4 deuterons in barrier penetration and fusion process, 
so that 4d simultaneous fusion should take place predominantly. The portion of 2D 
(usual) fusion rate is considered to be negligible.” 



1 



In this respect it can be noted that there are recent reports suggesting that hydrogen cluster can 
get reaction at very low temperature, forming the condition of superfluidity [5]. This seems to 
happen too in the context of Takahashi TSC condensate dynamics. Other study worth mentioning 
here is one that discussed molecular chessboard dynamics of deuterium [6]. 

The difference between this proposed study and recent work of Takahashi based on Langevin 
equation for cluster deuterium is that we focus on solution of Schrodinger-Langevin equation 
[7] [8] with PT-Symmetric periodic potential as we discussed in the preceding paper and its 
Gamow integral. The particular implications of this study to deuteron cluster will be discussed 
later. 

Another differing part from the previous study is that in this study we will also seek clues on 
possibility to consider this low probability problem as an example of self-organized criticality 
phenomena. In other words, the time required before CMNS process can be observed is actually 
the time required to trigger the critical phenomena. To our present knowledge, this kind of 
approach has never been studied before, although self-organized criticality related to Schrodinger 
equation approximation to Burger’s turbulence has been discussed in Boldyrev [12]. Nonetheless 
there is recent study suggesting link between diffusion process and the self-organized criticality 
phenomena. 

The result of this study will be useful to better understanding of anomalous phenomena behind 
Condensed matter nuclear science. 



Schrodinger-Langevin equation 

The Langevin equation is considered as equivalent and therefore has often been used to solve the 
time-independent Schrodinger, in particular to study molecular dynamics. 

Here we only cite the known Langevin equation [7, p. 29]: 



clX t = p r dt 


(1) 


dp = -d x A 0 ( X,)dt + Kp t dt + dW, 4 tTK 


(2) 



Takahashi & Yabuuchi also used quite similar form of the stochastic non-linear Langevin 
equation [8] in order to study the dynamics of TSC condensate motion. 
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Schrodinger equation with PT-Symmetric periodic potential 

Consider a PT-Symmetric potential of the form [9] [10]: 

V = k 1 .svn(b.r ) , 

where 

m 

"J- i - 1 



( 3 ) 



(4) 



Hence, the respective Schrodinger equation with this potential can be written as follows: 
¥"(r) = -fc 2 (r).¥(r) (5) 



where 

9 ill 9 yvj 

k(r ) = —[E-V(r)] = —[E - k v sm(b.r)] (6) 

tr h 

For the purpose of finding Gamow function, in area near x=a we can choose linear 
approximation for Coulomb potential, such that: 

V(x) - E = -a(x - a), (7) 



Substitution to Schrodinger equation yields: 



T „ + 2ma T = o 

n 2 



which can be solved by virtue of Airy function. 



( 8 ) 



Gamow integral 

In principle, the Gamow function can be derived as follows [11]: 

^i + P(x)y = 0 (9) 

clx 

Separating the variables and integrating, yields: 

J— = J -P(x).dx (10) 
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or 



y.dy = exp(-J P(x).dx ) +C) 



( 11 ) 



To find solution of Gamow function, therefore the integral below must be evaluated: 




( 12 ) 



For the purpose of analysis we use the same data from Takahashi’s EQPET model [3], [4], i.e. 
b=5.6fm, and ro=5fm. Here we assume that E=Vb=0.257MeV. Therefore the integral becomes: 



By setting boundary condition (either one or more of these conditions): 

(a) at r=0 then V G =-V b — 0.257 MeV 

(b) at r=5.6fm then Vl=k, sin(Zrr) -0.257 =0.257Mev, therefore one can find estimate of m. 

(c) Using this procedure solution of the equation (11) can be found. 

The interpretation of this Gamow function is the tunneling rate of the fusion reaction of cluster of 
deuterium (for the given data) corresponding to Takahashi data [4], with the difference that here 
we consider a PT- symmetric periodic potential. 

The numerical study will be performed with standard package like Maxima etc. Some plausible 
implications in Cosmology modeling will also be discussed. 
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